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^ . ABSTRACT 

o 

■ Context. The eclipsing low-mass X-ray binary XI 822-371 is the prototype of the accretion disc corona (ADC) sources. Its inclination 
angle 82.5°) is high enough that flux from the neutron star is blocked by the edge-on accretion disc. Because the neutron star's 
direct emission is hidden, its ADC emission is visible. 

Aims. We analyse two Chandra observations and one XMM-Newton observation to study the discrete features in this source and their 
variation as a function of the orbital phase, deriving constraints on the temperature, density, and location of the plasma responsible 
for emission lines. 

Methods. The HETGS and XMM/Epic-pn observed X1822-371 for 140 and 50 ks, respectively. We extracted an averaged spectrum 
and five spectra from five selected orbital-phase intervals that are 0.04-0.25, 0.25-0.50, 0.50-0.75, 0.75-0.95, and, finally, 0.95-1.04; 

■ the orbital phase zero corresponds to the eclipse time. All spectra cover the energy band between 0.35 and 12 keV. 
HH ' Results. We confirm the presence of local neutral matter that partially covers the X-ray emitting region; the equivalent hydrogen 
^; , column is 5 x 10"^ cm"^ and the covered fraction is about 60-65%. We identify emission lines from O vii, O viii, Neix Nex, Mgxi, 

(~ l . Mgxii, Si XIII, Si XIV, Fexxv, Fexxvi, and a prominent fluorescence iron line associated with a blending of Fei-Fexv resonant tran- 

^ [ sitions. The transitions of He-like ions show that the intercombination dominates over the forbidden and resonance lines. The line 

I ■ fluxes are the highest during the orbital phases between 0.04 and 0.75. 

O , Conclusions. We discuss the presence of an extended, optically thin corona with optical depth of about 0.01 that scatters the X-ray 

photons from the innermost region into the line of sight. The photoionised plasma producing the O viii, Ne ix, Ne x, Mg xi, Mg xii, 
^ ■ Sixm, and Sixiv lines is placed in the bulge at the outer radius of the disc distant from the central source of 6 x 10'° cm. The O vii 

and the fluorescence iron line are probably produced in the photoionised surface of the disc at inner radii. Finally, we suggest that the 
observed local neutral matter is the matter transferred by the companion star that was expelled from the system by the X-ray radiation 
pressure, which in turn originated in the accretion process onto the neutron star. (Abridged) 

> 

00' , 

■ 1. Introduction tween luminosity and spin-up rate, found that if X1822-371 is 
— »^ . at a luminosity of 10^*' erg s ' the neutron star magnetic field B 
O . X1822-371 IS a compact binary system with a period of 5.57 ^j^^^^j^ ^^^^ ^„ ^^jij^^jy ^^^^^^^^ of ~ 8 x 10'" G, whil e for a 

■ hr. The light curve of the source shows an almost sinusoidal luminosity of lO^*^ erg s ' B ~ 8 x lO'" G. iJonker et all (|2051 
^ ; modulation and partial eclipse. The partial echpse indicates a constrained the inclination angle of the source at 82.5° + 1.5° 

. highmchnation angle of the system. Initially, [Mason & Cordova| ^^j^g ultraviolet and visual data of X1822-371. This value im- 
> ■ (II982I) derived an inclmation angle between 76° and 84° fitting pij^^ ^^^^ modulation in the light curve of X1822-371 is not 
K/i : the light curve of X1822-371 in the infrared band. The mod- ^^^^g^ ^ip^ 1,^^^^^^ t^g^g can be ob served for inclination an- 
^ . ulation of the light curve is generally explained with the pres- gj^^ between 60° and 80° ffrarik et al.ll2002h . 

■ ence of a geometrically thick disc whose height varies depend- 

- - ■ ing on the azimuthal angle a nd occults part of the X-ray emission iBurderi et all (1201 Ol ) analysed the eclipse arrival times 
(White & Holt 1982; Mason & Cordovalll982l:lHellier & MasonI of X1822-371 using data from RXTE, XMM-Newton, and 
ri989l:iBavless et al. 2010). Chandra observations, spanning the years from 1996 to 2008. 

Analysing RXTE data of XI 822-37 l,|Jonker & van der KlisI Combining these eclipse arrival time measurements with those 
(1200 Ih detected a coherent pulsation at 0.59 s, and derived that already available ( covering the period from 1977 to 1996; see 
the orbit of the system is almost circular with an eccentricity iParmar etani2000h . the authors were able to tightly constrain 
less than 0.03, amass function of (2.03 +0.03) X 10"^ Mq, and a the orbital period derivative of the binary system to Porb = 
pulse period derivative of (-2.85 ± 0.04) x 10"'^ s s"' , indicating 1 .50(7) x 10"'° s/s, that is three orders of magnitude larger than 
that the neutron star is spiniiing up . The distance to XI 822-371 is expected from conservative mass tr ansfer driven by ma gnetic 
2.5 kpc dMason & Cordoval 19821) and the obse rved lumin osity i s braking and gravitational radiation. iBurderi etakl (1201 01) con- 
10^" erg s ' rsee e.g. 'Mason & Cordova'l982V laria et al.ll200lh . eluded that the mass transfer rate from the companion star is 
However, Jonker & van der Klis (2001), using the relation be- between 3.5 and 7.5 times the Eddington limit (i. e. 1 .8 x 10"^ 

Mo yr"' for a neutron star of 1.4 M©), suggesting that the mass 

Send offprint requests to: R. laria, e-mail: iaria(afisica.unipa. it transfer has to be highly non-conservative, with the neutron star 
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accreting at the Eddington limit and the rest of the transferred 
mass expelled from the system. 

'Bavless et al. (2010), studying the optical and UV data of 
XI 822-371, derived the new optical ephemeris for the source, 
finding a rapid change of the orbital period with PjP - (3.0 + 
0.3) X 10*" yr, similar to the value of the recent X - ray ephemeris 
of X1822-371 obtained by iBurderi et af] (l2010h . iBavless et all 
iOlO) showed that the accretion rate onto the neutron star 
should be 6.4x 10"^ yr"' in a conservative mass transfer, sug- 
gesting again a highly non-conservative mass transfer. Finally 
they fitted the optical light curve of XI 822-371, concluding that 
the disc has a vertically extended component that is optically 
thick at the optical wavelengths with a height of 0.5 times the 
accretion disc radius. The authors identified this component as 
the base of the disc wind. 

Recently, laria et all (1201 lb extended the work of 
IBurderi etal J (l2010l) . including a larg er sample o f data, 
and combined the optical/UV data from iBavless et al" ( 20101) 
to ob tain an updated ephemeris of XI 822-371. laria et all 
(1201 Ih obtained Porb = 1.59(9) x lO^^" s/s from the com- 
bined X-ray/optical/UV data set, consistent with the previous 
measurements. A similar value of the orbital period derivative 
(Pmh = 1.3(3) X 10""^ s/s) was independently obtained by 
iJain et all (1201 Oh by analysing data in the X-ray band. 

Because of the high inclination angle of the system and the 
partial eclipse, it has been argued that we do not observe the di- 
rect X-ray emission produced in the inner region of the system, 
but the emission coming from an extended cor ona above the 
disc, the so-called accretion disc corona (ADC, IWhite & Holt! 
|1982|) . The properties of the ADC are still debated. In a scenario 
with an optically thin ADC we should expect that only a small 
fraction of the emission produced in the inner region is scattered 
by t he ADC into the line of sight (see e.g. ' McClintock et alj 
1198 2: Hellier & Mason 1989). On the other hand, an optically 
thick coro na would imply an isotropic re-emission of the pri- 
mary flux (IWhite&Holtlll982l) . 

The X-ray spectrum of XI 822-371 has been widely debated. 
iHellier & Mason* ('1989') fitted the 1-10 keV EXOSAT spectrum 
of XI 822-371 with a black body peaked at 2 keV plus a flat 
power-law component and an iron emission line at 6.7 keV. The 
authors inferred that the area of the emitting black-body surface 
was smaller (~ 1/400) than the neutron star surface. They pro- 
posed that the black-body emission originates from the neutron 
star surface but the high inclination angle of the system and the 
tall outer disc would block part of the direct emission from the 
neutron star, and only a small percentage of the radiation is vis- 
ible to the observer because it is scattered by an optically thin 
corona. The power law component was associated with the emis- 
sio n from the inner accretio n disc. 

iHeinz & NowakI (l200lh analysed simultaneous observations 
with Rossi X-Ray Timing Explorer (RXTE) and ASCA of 
XI 822-371. They showed that both the source spectrum and 
light curve can be well-fitted by two equivalent models, repre- 
senting the scenario in which XI 822-371 has an optically thick 
corona or an optically thin corona, respectively. In the first case, 
no soft thermal component from the inner region contributes to 
the source spectrum, and the emission from the optically thick 
corona is described by a cutoff power law partially absorbed by 
a cold gas that is the atmosphere of the outer disc. In the sec- 
ond case, the model consists of a black-body component emitted 
from the central region, the neutron star surface and/or the inner 
disc, plus a Comptonised component produced in the optically 
thin corona and fitted with a cutoff power law in the spectrum. 
Both these models were able to describe the data adequately. 



Ilaria et al] (l200lh analysed a BeppoSAX observation of 
XI 822-371, finding that the continuum emission is well-fitted 
using a Comptonised component with an electron temperature 
of 4.5 keV that is partially covered by a local neutral absorber 
The authors suggested that the Comptonised emission is caused 
by a diffuse emission from the ADC while the partial cover- 
ing is produced by a cold wall subtending an angle of 16° with 
respect to the equato rial plane and place d in the external edge 
of the accretion disc. iParmar et alj ( l2000l) fitted the BeppoSAX 
data of XI 822-371 adopting a continuum model composed of a 
Comptonised component plus a black body. The electron tem- 
perature of the Comptonising region was 4.53 ± 0.02 keV, the 
optical depth was 26.2 + 0.6, the seed-photon temperature was 
0.15+ 0.02 keV and flie black-body temperature was 1 .27 + 0.03 
keV^ 

ICottam et al.l (1200 lb . analysing a Chandra observation 
of XI 822-371 with the High-Energy Transmission Grating 
Spectrometer (HETGS), identified several emission lines pro- 
duced in a photoionised plasma. The lines were associated with 
highly ionised iron, silicon, magnesium, neon, and oxygen, and 
with a prominent fluorescence iron line at 6.4 keV. The emis- 
sion region of the photoionised lines was placed at the inner part 
of the X-ray illuminated bulge that is at the point of impact be- 
tween the disc and the accretion stream from the companion star, 
whereas the origin of the fluorescence iron line was located in an 
extended region on the disc illuminated by the light scattered by 
the corona. 

Recently iJi et al.l (1201 ll) . analysing Chandra data sets also 
used for this paper, measured the Doppler velocities of several 
lines. They divided the observation into 40 intervals and moved 
a phase window with a width of 0.16 in phase across the 40 in- 
tervals, so that each interval had a phase shift of 0.025 from the 
previou s, and therefore the intervals partially overlap in phase. 
IJi etal.l (l2011). measuring the Doppler velocity modulation of 
several ionised lines, inferred that the observed photo-ionised 
emission lines originate from a confined region in the outer edge 
of the accretion disc near the hot spot. The photo-ionised plasma 
where the lines originate was consistent with ionisation param- 
eters ^ > 100. Ji et al. (2011), combining the disc size and rea- 
sonable assumptions for the plasma density (between 10'^ and 
10''* cm"-'), suggested that the illuminating disc luminosities are 
more than an order of magnitude higher than what is actually 
observed and concluded that the central emitting X-ray source is 
not directly observed. 

We follow a different approach by obtaining information on 
the lines as a function of the orbital phase from five independent 
intervals. In this paper we analyse two Chandra/HETGS obser- 
vations and one XMM-Newton observation of XI 822-371 for a 
total observing time of 140 and 50 ks, respectively. We discuss 
the averaged spectrum and the phase-resolved spectra and com- 
pare our results with past results discussed in literature. 

2. Observations 

2.1. Chandra observations 

X1822-371 was observed with the Chandra observatory from 
2008 May 20 22:53:00 to 2008 May 21 17:06:51 UT (Obs. ID. 
9076) and from 2008 May 23 13:20:56 to 2008 May 24 12:41:54 
UT (Obs. ID. 9858) using the HETGS for a total observation 
time of 66 and 84 ks, respectively. Both observations were per- 
formed in timed faint mode. The HETGS consists of two types 
of transmission gratings, the medium-energy grating (MEG) and 
the high-energy grating (HEG). The HETGS provides high- 
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Fig. 1. First-order MEG+HEG light curves of XI 822-371. The hght curves correspond to Obs. ID. 9076 (left panel) and 9858 (right 
panel), respectively. The bin time is 100 s for both light curves. 
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Fig. 2. Folded light curves correspondin g to Obs. ID. 90 76 (left panel) and Obs. ID. 9858 (right panel). The data were folded 
adopting the X-ray ephemeris obtained bv Jaria et alj (1201 li) and using 128 channels per period. The vertical lines indicate the phase 
intervals used for the phase-dependent spectral analysis (see section [331 ). 



resolution spectroscopy from 1.2 to 31 A (0.4-10 keV) with a 
peak spectral resolution of A/AA ~ 1000 at 12 A for the HEG 
first order. The dispersed spectra were recorded with an array of 
six charge-coupled devices (CCDs) that are part of the advanced 
CCD imaging spectrometer-S (Garmire et al. 2003). The current 
relative accuracy of the overall wavelength calibration is 0.05%, 
leading to a worst-case uncertainty of 0.004 A in the first-order 
MEG and 0.006 A in the first-order HEG. The brightness of the 
source required additional efforts to mitigate "photon pile-up" 
effects of the zeroth-order events. A 512 row "subarray" (with 
the first row =1) was applied during the observation, reducing 
the CCD frame time to 1.7 s. Pile-up distorts the count spec- 
trum because detected events overlap and their deposited charges 



are collected into single, apparently more energetic, events. We 
omitted the zeroth-order events from our spectral analysis and 
studied only the grating first-order spectra, which were not af- 
fected by pile-up. 

We processed the event list using available software 
(FTOOLS ver. 6.2 and CIAO ver. 4.1.2 packages) and com- 
puted aspect-corrected exposure maps for each spectrum, which 
allowed us to take into account their effects on the effective 
area of the CCD spectrometer We barycentred the events us- 
ing the tool axbary in CIAO and show the 100 s bin time light 
curves (merged HEG and MEG first-order events) corresponding 
to Obs. ID. 9076 and 9858 in Figure [1] In Fig. Ill we show the 
folded the light curves adopting the recent X-ray ephemeris of 
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llaria etalJ(l201lh and using 128 channels per period. The mean 
count rate of both observations is 5 cts/s out of the eclipse. 



2.2. XMM-Newton observation 



Folded period: 20054.2056 s 



> 



The XMM-Newton Observatory dJansen et all 1200 ih includes 
three 1500 cm^ X-ray telescopes each with an European Photon 
Imaging Camera (Epic, 0.1-15 keV) at the fo cus. Two of 
the EPIC imaging spectrometer s use MPS CCDs (iTurner et al.l 
l200lh and one uses pn CCDs ( Struderet al.ll2001l). Reflection 
gratin g spectrometers (RGS, 0.35-2.5 keV, den Herder et al.l 
1200 ih are located behind two of the telescopes. The region of 
sky containing XI 822-371 was observed by XMM-Newton be- 
tween 2001 March 07 13:12:48 UT and March 08 03:32:53 UT 
(Obs. ID. 01 1 1230101) for a duration of 53.8 ks. During the ob- 
servation the MOSl and MOS2 camera were operated in fast 
uncompressed mode and small window mode, respectively. The 
Epic-pn camera was operated in timing mode with medium fil- 
ter during the observation. In this mode only one central CCD 
is read out with a time resolution of 0.03 ms. This provides a 
one dimensional (4'.4 wide) image of the source with the sec- 
ond spatial dimension being replaced by timing information. The 
faster CCD readout results in a much higher count rate capabil- 
ity of 800 cts/s before charge pile-up becomes a serious problem 
for point-like sources. The Epic-pn count rate of the source was 
about 50 cts/s, which avoids telemetry and pile-up problems. 

For our analysis we used only RGS and Epic-pn data since 
the Epic-pn covers the same energy range as the MOS CCDs and 
has a larger effective area. We extracted the X-ray data products 
of RGS and Epic-pn camera using the science analysis software 
(SAS) version 9.0.0 extracting only single and double events 
(patterns to 4) from Epic-pn data. We selected source Epic- 
pn events from a 69.7" wide column (RAWX between 29 and 
45) centred on the source position (RAWX= 38). 

Using the SAS tools, we accumulated the RGS 1 and Epic- 
pn background light curves. The RGSl background light curve, 
plotted in Fig. [3] (left panel), shows a flare (possibly due to so- 
lar contamination) during the first 10 ks. We excluded this time 
interval from our analysis because of its spurious origin. The 
Epic-pn background light curve was extracted from the columns 
between 2 and 18 and for energies between 10 and 12 keV (see 
Fig. [3] right panel). Together with the large flare in the first 10 ks, 
a smaller flare is seen between 37 and 41 ks in the Epic-pn light 
curve; we excluded these two time intervals and barycentred the 
events using the tool barycen in SAS to study the orbital-phase 
resolved spectrum. 

Finally, we show in Fig. |4]the Epic-pn folded light curves 
selected in the energy bands 1-3 keV and 3-10 keV, and the 
folded hardness ratio (HR) obtained from the ratio of the 3-10 
keV folded light curve over the 1-3 keV fold ed hght curve. We 
note that the HR is similar to that obtained bv ' laria et al.l (1200 ih 
and Parmaretal. (2000) using BeppoSAX data of XI 822-371. 
The largest effective area of the Epic-pn allows us to clearly ob- 
serve a small variation of the HR, from 1 .5 to 1 .9, which imme- 
diately preceeds the eclipse at the orbital pha se 0.949. This fea - 
ture was akeady present in the HR shown bv laria e t al.l (1200 ih . 
[Parm ar et al. (2000), and Heflier & Mason ( 1989). Furtiiermore, 
we confirm that the depth of the partial eclipse in creases with 
increasing energy as shown bv lParmar et al.l (l2000l) . 
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Fig. 4. Upper (middle) panel: Epic-pn 1-3 keV (3-10 keV) 
folded light curves. Lower panel: Epic-pn hardness ratio for 
the two energy-selected light curves shown above. The data 
were folded adop ting the recent X-ray ephemeris obtained by 
llaria et al] (1201 lb and using 256 channels per period. 



3. Spectral analysis 

3.1. The averaged spectrum 

Chandra data were extracted from regions around the grating 
arms; to avoid overlapping between HEG and MEG data, we 
used a region size of 25 and 33 pixels for the HEG and MEG, re- 
spectively, along the cross-dispersion direction. The background 
spectra were computed, as usual, by extracting data above and 
below the dispersed flux. The contribution from the background 
is 0.4% of the total count rate. We used the stan dard ClAO tools 
to create detector response files ( lDavisll2001h for the HEG-1 
(MEG-1) and HEG+1 (MEG -hi) order (background-subtracted) 
spectra. After verifying that the negative and positive orders 
were compatible with each other in the whole energy range, 
we coadded them using the script add^ratingjorders in the 
ClAO software, obtaining the first-order MEG spectrum and 
the first-order HEG spectrum for each of the two observations. 
We summed the first-order HEG (and MEG) spectra of the two 
Chandra observations to increase the statistics with the ClAO 
script add-gratingjpectra and obtained a total exposure time of 
142 ks. The MEG and HEG energy range is 0.5-6 keV and 0.8- 
10 keV, respectively. Finally, we rebinned the first-order MEG 
and first-order HEG spectra to have at least 25 counts per energy 
channel. 

We extracted the RGSl, RGS2, and Epic-pn spectra from 
the whole observation with an exposure time of 43, 42, 38 ks, 
respectively. We adopted a 0.35-2 keV energy range for RGS 
spectra and 0.6-12 keV for the Epic-pn spectrum. The Epic-pn 
spectrum was rebinned so that the energy resolution was not 
oversampled by more than a factor 4 and to have at least 25 
counts per energy channel. The RGSl and RGS2 spectra were 
rebinned to have at least 25 counts per energy channel. 

The HETGS and XMM spectra were fitted separately us- 
ing XSPEC version 12.6.0. We fitted the co ntinuum emis- 
sion adopting the model proposed by llaria et al. (2001) for 
the broad band BeppoSAX spectrum of X1822-371: i.e. a 
Comptonised component (CompTT in XSPEC) partially ab- 
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Fig. 3. RGS 1 background light curve (left panel) and Epic-pn background light curve (right panel). The horizontal segments indicate 
the time interval excluded from our analysis. The bin time is 100 s. 



Table 2. Averaged spectrum. Emission lines of XI 822-371. 





Identification 


Expected 


E 


cr 


I 


eq. width 


significance 






(keV) 


(keV) 


(eV) 


(10^^ ph cm-- s-') 


(eV) 


(0-) 




Ovii(i) 


0.5687 


0.5684 + 0.0003 


1.1 ±0.2 


34 ±6 


19 ±6 


6 




O VIII 


0.6536 


6529+''-'"'2 


1.1 ±0.3 


9±2 


6±2 


5 




Ne IX (i) 


0.9149 


0.9147 ±0.0003 


1.7 ±0.3 


4.3 ±0.5 


5.2 ± 1.0 


9 




Nex 


1.0218 


1.0212 ±0.0003 


1 74+0.14 

^ - ' -0.30 


2.4 ± 0.3 


3.4 ±0.8 


8 


HETGS 


Mgxi (i) 


1.3433 


1.3430 ±0.0004 


2.8 ±0.5 


1.40 ±0.14 


2.8 ±0.6 


10 


Mgxii 


1.4723 


1.4706 ±0.0005 


2.5 ± 0.4 


0.82 ±0.10 


1.9 ±0.6 


8 




Si xm (i) 


1.8542 


1.8551 ±0.0010 


3 (fixed) 


0.46 ± 0.08 


1.4 ±0.6 


6 




Si XIV 


2.005 


2.0043 ± 0.0007 


3.7 ±0.6 


0.74 ± 0.09 


2.5 ±0.7 


8 




Fei 


6.404 


6.393 ± 0.002 


18 ±3 


3.2 ±0.2 


40 ±6 


16 




Fe xxv(i) 


6.675 


6.675 (fixed) 


13 (fixed) 


0.6 ± 0.2 


8±4 


3 




Fe XXVI 


6.966 


6.952 ± 0.004 


13 ±6 


1.5 ±0.2 


22 ±8 


8 


RGS 


Ovii(i) 
vm 
Ne IX (i) 


0.5687 
0.6536 
0.9149 


56869+"-™" 
651794+" ™"'^ 
0.9139 ±0.0011 


1.4 ±0.2 

<0.7 
4.4 ± 1.3 


23.2 ± 1.3 
3.6 ± 0.4 
3.8 ±0.3 


22 ±2 
4.6 ± 1.0 
8.0 ± 1.3 


18 
9 
13 




Nex 


1.0218 


1.0218 (fixed) 


2 (fixed) 


1.3 ±0.2 


3.2 ± 1.1 


7 




Mgxi(i) 


1.3433 


1.3433 (fixed) 


3 (fixed) 


0.6 ± 0.2 


2.0 ± 1.1 


3 




Si XIV 


2.005 


1 084+0.010 


4 (fixed) 


1.49 ±0.15 


8±2 


10 


Epic-pn 


Fei 
Fexxv (i) 


6.404 
6.675 


6.402 ± 0.004 
6.703 ±0.013 


75 ±8 
20 (fixed) 


3.24 ±0.13 
0.70 ± 0.08 


56 ±5 
12 ±3 


25 
9 




Fe XXVI 


6.966 


6.991 ±0.009 


56+** 


1 oo+O.U 


26 ±5 


23 



Note — Best-fit values of the HETGS and XMM emission lines in X1822-371. Uncertainties are at the 68% c. 1. for each parameter. In the 
first column we report the atomic transition, from the second to the fifth column the expected rest-frame energy, the centroid, the width, and the 
intensity of each line. In the sixth and seventh column we report the equivalent width and the significances, in units of cr, of each line. 



sorbed by neutral matter and absorbed by in terstellar mat- 
ter (phabs in XSPEC, using the abundances of lAsplund et al.l 
(^09)). We also took into account the effect of Thomson 
scattering in the local cold absorber using the cabs compo- 
nent in XSPEC. The adopted model consists of the component 
phabs*(f"cabs"phabs''CCompTT)-i-Cl-£)"CompTT), where f 
is the neutral matter covering fraction; the first phabs compo- 
nent was used to fit the photoelectric absorption by interstellar 
neutral matter, the second phabs component takes into account 



the photoelectric absorption by local neutral matter. We imposed 
that the value of the equivalent hydrogen column density of the 
cabs component is the same as that of the neutral local matter, 
assuming that the local neutral matter is responsible for both the 
photoelectric absorption and Thomson scattering. 

Fitting the RGS and Epic-pn spectra, we obtain ax^^J^d-o.f.) 
of 1.86(2422). The residuals reported in Fig.|5](the second panel 
from the top) show several local features and a mismatch be- 
tween the RGS and Epic-pn data in the energy range between 0.6 
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Table 1. Averaged spectrum. Best-fit values of the continuum 
parameters. 



Parameters 

N„ (10-^ cm-^) 

N„,,^, (10^2 ^^-2) 

/ 

kTifc (keV) 
Nfci (xlO-3) 

kTo (keV) 
kT, (keV) 

T 

^CompTT (XlQ--) 



HETGS 



XMM 



0.163 (fixed) 0635™"" 

4.84 ±0.10 4.97 ±0.12 
0.654 ± 0.012 0.607 ±0.009 



0.061 (fixed) 
1.5 ± 1.0 

0.061 (fixed) 

■^•"'-0.14 

19.1 ±0.7 
12.7 ± 0.4 

0.82(6125) 



0.061 ± 0.003 
1.0 ±0.2 

0.061 ± 0.003 
3.05 ± 0.04 
21.0 ±0.4 
7.7 ±0.2 

1.06(2400) 



Note — Uncertainties are at the 90% c. 1. for each parameter. The pa- 
rameters are defined as in XSPEC. The model is discussed in the text. 
The equivalent hydrogen column of the local neutral matter is indicated 
with ^Hpr' 'he covered fraction of the emitting surface is indicated with 
/. The. xled^d.o.f.) values are obtained taking into account the emission 



lines reported in Table |2] 




1 2 

Energy (keV) 

Fig. 5. Data and residuals of RGSs and Epic-pn spec- 
tra. Second panel from the top: residuals using the model 
phabs* Cf *cabs*phabs " (CompTT) -i-C l-£) -CompTT) . Third 
panel from the top: residuals adding nine Gaussian components 
to the previous model to fit the lines. Bottom panel: residuals 
adding a black-body component to the model. 



and 0.8 keV. We improved the model by adding nine Gaussian 
components to fit the localised features. The centroids of the 
Gaussian profiles are at 0.569, 0.652, 0.914, 1.02, 1.34, 1.98, 
6.40, 6.70, and 6.99 keV, associated with the Ovii intercombi- 
nation line, O viii, Ne ix intercombination line, Ne x, Mg xi inter- 
combination line. Si xiv, Fe i, Fe xxv intercombination line, and 
Fexxvi, respectively. The addition of the Gaussians improved 
the fit, giving a A;^' of 1691 and a xl^iiA-O-^) of 1.17(2401); 




channel energy (kcV) 



liiiitil 



0.96 0.98 
channel energy (keV) 




channel cnci'iiy (kcV) 



Fig. 6. Residuals of first-order MEG (black points) and HEG 
(red points) spectra using the continuum shown in the text. The 
emission lines are indicated in the panels. 
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the corresponding residuals are shown in Fig.|5](the third panel 
from the top). To improve the fit at low energies, where RGS 
and Epic-pn data are mismatched, we added a black-body com- 
ponent (bbody in XSPEC) to our model, partially absorbed like 
the Comptonised component. Since the best-fit gives a value of 
0.06 keV for the black-body temperature and for the seed-photon 
temperature of the CompTT component, we constrained the two 
parameters to assume the same value and fitted the data again. 
The addition of the black-body component improved the fit, giv- 
ing a l^x^ of 250 and a ;|^2^j(d.o.f.) of 1.06(2400). We report the 
corresponding residuals in Fig.|5](bottom panel). The best-fitting 
values of the parameters of the continuum are reported with an 
associated error at 90% confidence level (c.l.) in Table[T|and the 
best-fitting values of the line parameters are reported in Table |2] 
(associated error at 68% c. 1.). 

We fitted the HETGS data with the same best-fit continuum 
as adopted above, obtaining a value of;t'^^^(d.o.f.) = 0.95(6155); 
we observed several emission features in the residuals as shown 
in Fig.|6l To fit the residuals, we added eleven Gaussian compo- 
nents and identified the emission lines associated with the fol- 
lowing transitions: O vii intercombination line, O viii Lyo', Ne ix 
intercombination line, Ne x Lyor, Mg xi intercombination line, 
Mg xn Lyo', Si xiii intercombination line. Si xiv Lya, Fe i Ka, 
Fe xxv intercombination line, and, finally, Fe xxvi Lya. The ad- 
dition of the Gaussian components improved the fit with a t^x^ 
of 807 and a value of ^(d.o.f .) of 0.82(6125). The best-fit val- 
ues of the continuum parameters and the best-fit values of the 
line parameters are reported in Tables [T]and|2] respectively. 

The black-body component added to fit the XM M spectrum 
below 0.8 keV is not present in the model shown bv llaria et alJ 
(1200 lb . This component was also added to fit the Chandra av- 
eraged spectrum, fixing its temperature to the best-fit value ob- 
tained by the XMM averaged spectrum. We obtained a consistent 
value for its normalisation. 

Finally, we checked the presence in the RGS spectrum of the 
Mgxn and Sixiii lines, which are significantly detected in the 
HETGS spectrum. We added two Gaussian components in the 
XMM averaged spectrum, keeping the centroids and the widths 
fixed to the values obtained from Chandra. We obtained an up- 
per limit to the flux of the Mg xii and Si xiii line of < 0. 12 x 10""^ 
and < 0.5 x 10"^ ph cm"^ s"', respectively. These fluxes suggest 
that the Mg xii line might be absent during the XMM observa- 
tion while the Si xiii line is compatible with that observed in the 
Chandra averaged spectrum. 

3.2. Helium-like series and the Ne ix RRC feature 

The intercombination lines associated with Ovii (RGSl spec- 
trum), Ne IX (MEG spectrum), and Mg xi (MEG spectrum) have 
a significance larger than 7 cr but we do not have evidence for 
the corresponding forbidden and resonance lines. To determine 
the intensity of the O vn lines, we analysed the energy range be- 
tween 0.555 and 0.58 keV, fixing the best-fitting values of the 
continuum emission. We studied the O vii triplet using RGS 1 and 
MEG data. In the MEG and RGS 1 spectra we fixed the energy of 
the forbidden line at the rest-frame value of 0.5610 keV, while 
the energies of the intercombination and resonance lines were 
let free to vary. The line widths were constrained to assume the 
same value. Using the best-fitting intensities reported in Table 
[3] we derived the parameter G - i f + i)/r and R - f/i, where 
/, /, and r are the intensities of the forbidden, intercombination, 
and resonance lines, respectively. We find that G = 4.3^j g and 
G = 7 + 4 for MEG and RGSl data and R is less than 0.07 and 
0.02 for MEG and RGSl, respectively (the associated errors are 




0.555 0.56 0.565 0.57 0.575 0.58 

channel energy (keV) 




0.555 0.56 0.565 0.57 0.575 0.58 

channel energy (keV) 













^ 1 1 1 1 1 1 _ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 





0.905 0.91 0.915 

channel energy (keV) 




1.33 1.34 1.35 1.36 1.37 

channel energy (keV) 

Fig. 7. Forbidden (green), intercombination (blue), and reso- 
nance (red) lines of the O vii, Ne ix and Mg xi triplets. From top 
to bottom: RGS 1 data and residuals of the O vii triplet; first-order 
MEG data and residuals of the O vii triplet; first-order MEG data 
and residuals of the Ne ix triplet; first-order MEG data and resid- 
uals of the Mg XI triplet. 



at 68% c. 1.). We used the MEG data between 0.89 and 0.93 keV 
to estimate G and R of the Ne ix triplet. We fixed the forbidden 
line energy at the rest-frame value of 0.9051 keV and imposed 
the widths of the three lines to be the same, finding G 



4 4+3.9 



and R = 0.04+° °^. FinaUy, we used the MEG data between 1.32 
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and 1 .37 keV to determine G and R of the Mg xi triplet, we fixed 
the energies of the forbidden and resonance line at the rest-frame 
values of 1.3312 and 1 .3523 keV, while the line widths were con- 
strained to assume the same value. We obtained G - 3.7;^j j and 
R ' 



ll+o ''^ 



Table 3. Best-fit values of the O vii, Ne ix, and Mg xi triplets. 





E cr 
(keV) (eV) 


I 

(lO^" phcm-' s-') 




Ovii 




f 


0.5610 (fixed) 1.0 ±0.2 


< 2 


MEG i 


0.5684 ± 0.0002 1.0 ±0.2 


26 ±3 


r 


0.5738 ± 0.0010 1.0 ±0.2 


6±2 



RGSl 



MEG 



MEG 



0.5610 (fixed) 
0.5725+f ■ 



1.0: 
1.0: 
1.0: 

Neix 



0.9051 (fixed) 
0.9145 ±0.0002 
0.9194 ±0.0008 

Mg 

1.3312 (fixed) 
1.3429 ±0.0004 
1.3523 (fixed) 



1.3: 
1.3: 
1.3: 
XI 

2.5 : 
2.5: 
2.5: 



0.2 
0.2 
0.2 

0.2 
0.2 
0.2 

0.4 
0.4 
0.4 



< 0.2 

12.6 ± 1.3 

I 9+1.0 



12*°-^^ 

"•^ -0.12 



3.4: 
0.8: 



0.4 
0.3 



12+°'2 
1.10 ±0.14 

33+° " 



Note — Centroid, width, and intensity of the He-like transition lines. 
The errors associated with each parameter are at 68% c. 1.. 



ICottam et al.l 12001) found a radiative recombination contin- 
uum (RRC) feature associated with the Ne ix transition in a pre- 
vious Chandra/HETGS observation of XI 822-371. They found 
an energy threshold of 13 ± 7 eV, corresponding to a plasma 
temperature of (1.5 + 0.8) x 10^ K. We searched for this fea- 
ture in our data by adding an inverse edge component (redge 
in XSPEC) to assess the presence of the RRC feature. We fixed 
the energy threshold at 1 . 196 keV (that is the energy of the Ne ix 
transition) and found only marginal evidence for the RRC fea- 
ture {Ax^ - 2.8), with a temperature <14 eV corresponding to a 
temperature of < 1.6x1 0^ K. 

3.3. Phase-dependent spectral analysis 

We selected the following five phase-intervals 0.04-0.25, 0.25- 
0.50, 0.50-0.75, 0.75-0.95, and 0.95-1.04, respectively. This se- 
lection allows us to constrain both continuum emission and dis- 
crete features for each of phase-interval. The intervals corre- 
spond to the count rate rise after the eclipse, the maximum, the 
count rate decrease, the minimum and, finally, the eclipse, re- 
spectively, as shown in Fig.|2] 

For each phase-interval we produced a spectrum using the 
same procedure as adopted for the averaged spectrum. We label 
the spectra extracted from the phase-intervals 0.04-0.25, 0.25- 
0.50, 0.50-0.75, 0.75-0.95, and 0.95-1.04 with 1, 2, 3, 4, and 5, 
respectively. The exposure times of the spectra corresponding 
to the five phase-intervals are reported in Table |4] The MEG, 
HEG, RGS, and EPIC-pn selected energy ranges are the same as 
adopted to fit the averaged spectrum. The model to fit the XMM 
spectra is the same as used to fit the averaged spectrum, while we 
fitted the five Chandra spectra excluding the black-body compo- 
nent because of the low statistics below 0.6 keV in the MEG. 
Furthermore, we fixed the Nh value to 0.163 x 10^^ cm"^ in the 



> 



0.04-0.25 



0.25-0.50 0.50-0.75 
Orbital Phase 



0.75-0.95 



0.95-1.04 



Fig. 8. Parameters of the Comptonised component obtained fit- 
ting the XMM data: electron temperature (top), optical depth 
(middle), and the unabsorbed flux in units of 10"'° erg s"' cm"^ 
(bottom) in the 0.4-12 keV energy band, the open squares and 
the dots are the values obtained from the Chandra and XMM 
spectra, respectively. 



Chandra and XMM spectra and, in the Chandra spectra, the seed 
photon temperature of the Comptonised component was fixed to 
the values obtained from the best-fit of the corresponding XMM 
spectrum. 

We obtained acceptable fits, the XMM spectra having 
xl/d.o.f.) of 1.00(815), 1.04(982), 0.98(1080), 1.00(707), and 
0.97(475) for spectra 1, 2, 3, 4, and 5, respectively. The best-fit 
parameters are shown in Tables |5] and |6l The electron tempera- 
ture kTe does not change with the orbital phase (see top panel 
of Fig. |8]l, while the optical depth t of the Comptonised com- 
ponent is about 20 at all orbital phases except between 0.75- 
0.95, when it increases up to 23 (see middle panel of Fig. [8). 
We show the unabsorbed flux in the 0.4-12 keV energy range 
associated with the Comptonised component that was obtained 
using the cflux component of XSPEC in the bottom panel of 
Fig. [8] where the open squares and the dots indicate the Chandra 
and XMM unabsorbed fluxes, respectively. The behaviour of the 
unabsorbed fluxes obtained from the XMM and Chandra phase- 
resolved spectra is similar along the orbit. The Chandra unab- 
sorbed fluxes are a factor 1 .4 higher than the those of XMM in 
all selected phase-intervals. 

The equivalent hydrogen column associated with the local 
neutral matter, Nh^,^., reaches the maximum at phases between 
0.25-0.50 when the unabsorbed flux is maximum. It assumes the 
minimum value at phases between 0.75-0.95 (see Fig |9l upper 
panel). A similar behav iour of A^h,,, and / was akeady found by 
iHeinz & NowakI (l200lh out of eclipse from analysing ASCA and 
RXTE data simultaneously (see Table 2 in their work). However, 
we found that the values of Nh and / increase again during the 
eclipse, while iHeinz & Nowakf(l200lf) found a further drop. 

The emission lines are clearly detected during most orbital 
phases, as shown in Table |6l The HETGS lines are evident dur- 
ing the phases between 0.04 and 0.75, that is out of the minimum 
in the light curve and eclipse. In spectra 1, 2, and 3 we fixed only 
the width of the O vii line at the values obtained in the averaged 
spectrum because of the low statistics below 0.6 keV. In spectra 2 
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Table 4. Exposure times of the selected phase-intervals. 



Phase interval 


Obs. 9076 


Obs. 9858 


RGSl 


RGS2 


Epic-pn 


0.04-0.25 


15.0 ks 


17.2 ks 


8.1 ks 


7.9 ks 


7.3 ks 


0.25-0.50 


14.7 ks 


20.8 ks 


9.8 ks 


9.5 ks 


13 ks 


0.50-0.75 


14.7 ks 


19.6 ks 


12.8 ks 


12.4 ks 


8.5 ks 


0.75-0.95 


11.7 ks 


15.7 ks 


7.9 ks 


7.6 ks 


7.2 ks 


0.95-1.04 


6.1 ks 


7.0 ks 


3.6 ks 


3.6 ks 


1.8 ks 



Note — We show the selected phase-intervals in the first column, the corresponding exposure times of the two HETGS spectra in the second and 
third columns, the exposure times of the RGSl, RGS2, and Epic-pn spectra in the fourth, fifth, and sixth columns, respectively. 

Table 5. Best-fit values of the continuum parameters as a function of the orbital phase of XI 822-371 



Spectra 


Parameters 


0.04-0.25 
1 


0.25-0.50 

2 


Phase intervals 
0.50-0.75 

3 


0.75-0.95 
4 


0.95-1.04 
5 




Nh (10^2 cm--) 


0.163 (fixed) 


0.163 (fixed) 


0.163 (fixed) 


0.163 (fixed) 


0.163 (fixed) 


HETGS 


N„„^ (10^2 cm-2) 
/ 


4.9 ± 0.2 


4.8 ±0.2 
0.63 ± 0.03 


4.7 ±0.3 
0.65 ± 0.04 


4.8 ±0.3 
0.66 ± 0.04 


4 6+"-^ 

^•"-0.6 

65+" "^ 

"•"-^-0.09 


kTo (keV) 
kT, (keV) 

T 

Nc,„„,.„ (10-2) 


0.060 (fixed) 

T T 1+0.31 

-'■-'^-0.07 

17.80 + 0.14 

13 6+"^ 


0.059 (fixed) 

-0.3 

20 ±2 
17.8 ±0.8 


0.054 (fixed) 
2.8 ±0.3 
20 ±2 
13.2 ±0.9 


0.054 (fixed) 

T Q + 1.0 

-'•■'-04 

18 ± 2 
10.7:|;° 


0.059 (fixed) 

T 4+71.0 

18+t, 

-12 

8.6!",^ 




Unabs. Flux 




10.85!":?^ 


o Q7+0.05 
"•^'-0.15 


7 63+" "'' 

'■"-'-0.11 


6 32+" "' 

"■■-'^-0.18 




xL(d.o.f.) 


0.70(3687) 


0.73(4035) 


0.67(3611) 


0.63(3029) 


0.49(1467) 




N„ (10^2 cm--) 


0.163 (fixed) 


0.163 (fixed) 


0.163 (fixed) 


0.163 (fixed) 


0.163 (fixed) 


XMM 


N„^^ (10^2 cm-2) 
/ 

kT,* (keV) 
Ntt (xlO-3) 


5.1 ±0.2 
0.604!«;«M 

060+" '"''' 

"•"""-0,007 

1 0+" 

^•"-0.3 


5.5 ±0.2 
0.612 ±0.014 

0.059 ± 0.006 

1 9+0.8 
^■^-0.4 


5.1 ±0.2 

643+" "" 

"■"^■'-0.017 

0.054 ± 0.005 
1.4 ±0.7 


4.4 ± 0.2 
0.58^^^ 

1 6+'-' 

i.U_Q7 


5.1 ±0.4 

63+" °^ 

"■"-'-0.05 

0.059 ± 0.015 

7+'-'' 

"•'-0.5 




kTo (keV) 
kT, (keV) 

T 

Nc<„„,„, (xlO-2) 


060+" 

"•"""-0.007 

3.06 ± 0.09 
20.3 ± 0.7 

7 7+0.5 
' • ' -0.3 


0.059 ± 0.006 
3.11 ±0.07 
20.6 ± 0.5 
8.8 ±0.5 


0.054 ± 0.005 
3.06 ± 0.07 
20.0 ± 0.6 
8.99 ±0.12 


054+" ""' 

"•"-"^-0.005 

2.98 ± 0.08 
22.5 ±0.8 
6.5 ± 0.3 


0.059 ±0.015 
3.0 ±0.2 
20.0+[-^ 
5.4 ± 0.7 




Unabs. Flux 


5.99 ± 0.03 


7.22 ± 0.03 


6.60 ± 0.03 


5.52 ± 0.03 


4 05+" "-"^ 






1.00(839) 


1.06(1015) 


0.97(1115) 


1.01(728) 


0.98(487) 



Note — Uncertainties are at the 90% c. 1. for each parameter. The 0.4-12 keV unabsorbed flux is given in units of 10 '" erg s ' cm ^. The 
xl^Jyd.o.f.) anA- x^{d-o.f.) values are obtained taking into account the emission lines shown in Tab.[6] 



and 3 we fixed the width of the O viii line at the best-fit value ob- 
tained in spectrum 1 . The other parameters of the O vii and O viii 
lines were left free to vary in spectra 1, 2, and 3. Furthermore, all 
parameters associated with the other detected lines were left free 
to vary in spectra 1, 2, and 3. In spectrum4 (phase 0.75-0.95) we 
fixed both the energies and the widths of the O vii and O viii by 
obtaining upper limits for the corresponding fluxes; moreover, 
we fixed the widths of the Ne x, Mg xi, Mg xii. Si xiv, and Fe xxvi 
lines at their best-fit values. In spectrum 5 (phase 0.95-1.04) we 
found an upper limit for the O vii and O vm, while we fixed the 
widths of the other detected lines at their best-fit values. 



In the XMM spectra we fixed the width of the weak Fe xxv 
line at 20 eV as in the averaged XMM spectrum. We fixed the 
width of the O vm line at the corresponding upper limit obtained 
from XMM spectrum 1 in spectra from 2 to 5. In XMM spectrum 
2 we fixed the width of the Fe xxvi lines at the corresponding 
upper limit obtained in spectrum 1 . In XMM spectrum 4 we fixed 
the energy and the width of the Ne ix line, the width of the Si xiv 
line, and the width of the Fe xxv line at their best-fit values. In 
XMM spectrum 5 the widths were fixed at their best-fit values. 
The equivalent widths in units of eV of the most prominent lines 
are shown in Table Q 
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Table 6. Emission lines of XI 822-371 





Id. 


E 




1(10-") 


significance 


E 


cr 


^lO-*) 


significance 






(keV) 


(eV) 


(ph cm"^ s"') 


(0-) 


(keV) 


(eV) 


(ph cm"^ s"') 


(0-) 








0.04-0.25 Phase-Interval 






0.25-0.50 Phase-Interval 






O VII (i) 
O VIII 


"■5690_o.ooi5 
0.6525 ± 0.0008 


1.1 (fixed) 
1.3 ±0.5 


57 ±23 
15 ±6 


2.5 
2.5 


n 5/:oQ+0.0029 

0.3t)SS„gQj, 

0.6538 ± 0.0012 


1.1 (fixed) 
1.3 (fixed) 


00+27 
-20 

10 ±5 


1.7 

2 




Ne IX (i) 
Nex 


0141+0.0002 
1.0207 ±0.0004 


0.9 ± 0.3 
1.2 ±0.4 


6.0 ± 1.2 
3.4 ±0.8 


5 

4.3 


0.9159 ±0.0005 
1.0216 ±0.0004 


1.8 ± 0.5 
1.7 ±0.4 


6.1 ± 1.3 
3.4 ±0.7 


4.7 
4.9 


HETGS 


Mg XI (i) 


1.3423 ±0.0004 


1.2 ±0.5 


1.5 ±0.3 


5 


1.3440 ±0.0006 


3 2+'-'' 
2 2+^'^ 


2.4 ± 0.4 


6 




Mgxii 


1.4705 ±0.0004 


1.3 ±0.5 


1.2 ±0.2 


6 


1.4725 ± 0.0009 


0.9 ± 0.3 


3 




Si XIV 


2.0038 ± 0.0008 


^■^-0.8 


1.1 ±0.2 


5.5 


2.0071 ±0.0011 


2.8 ± 1.0 


0.8 ± 0.2 


4 




Fe I 


6.396 ± 0.003 


14 ±3 


4.1 ±0.5 


8.2 


6.396 ± 0.004 


19 ±5 


4.4 ± 0.6 


7.3 




Fe XXVI 


6.960 ± 0.007 


< 22 


2.1 ±0.5 


4.2 


6 946+0 008 

°-^™-0.002 


< 9 


1.8 ±0.5 


3.6 




O VII (i) 
O VIII 


0.5684 ± 0.0003 
0.651 ±0.002 


1.6 ±0.5 

< 3 


25 ±3 
3.0 ± 1.0 


8.3 

3 


0.5688+U'iJilnT 

-0.0003 

0.6536 ± 0.0009 


1.1 ±0.3 
3 (fixed) 


28 ±3 
5.7 ± 1.2 


9.3 
4.8 


AMJVl 


Ne IX (i) 


0.913 ±0.002 


4±3 


4.3 ±0.7 


6.1 


0.9118 ±0.0014 


3 ±2 


5.6 ±0.6 


9.3 


Si XIV 


1.99 ±0.02 


44+20 


2.1 ±0.6 


3.5 


1.978 ±0.012 


56 ±20 


2.9 ±0.6 


4.8 




Fei 


6.403 ± 0.010 


64 ±25 


2.9 ±0.4 


7.3 


6.408 ± 0.007 


82± 10 


4.0 ± 0.3 


13 




Fe XXV (i) 


6.68 ± 0.03 


20 (fixed) 


1.0 ±0.2 


5 


6.75 ± 0.02 


20 (fixed) 


0.84 ±0.15 


5.6 




Fe XXVI 


6 99+0.02 


< 39 


1.5 ±0.2 


7.5 


7.005 ± 0.012 


39 (fixed) 


1.4 ±0.2 


7 








0.50-0.75 Phase-Interval 






0.75-0.95 Phase-Interval 






VII (i) 


5695+""'"" 


1.1 (fixed) 


48 ±20 


2.4 


0.569 (fixed) 


1.1 (fixed) 


< 77 






VIII 


6525+™'" 


1.3 (fixed) 


9±5 


1.8 


0.653 (fixed) 


1.3 (fixed) 


< 13 


- 




Ne IX (i) 


0.9142 ± 0.0005 


1.3 ±0.5 


3.7 ± 1.2 


3.1 


9138+" ™2 


2 4+2-' 


3±2 


1.5 


HETGS 


Nex 


1 0226+''"002 

1 .U^Z.U_Q QQQ2 


< 0.7 


1.4 ±0.5 


2.8 


1.0208 ± 0.0006 


0.7 (fixed) 


1.2 ±0.5 


2.4 




Mgxi (i) 


1.344 ±0.002 


5 l+^-O 

-"■'-1.4 


1.2 ±0.4 


3 


1.343 ±0.002 


3 (fixed) 


0.6 ± 0.3 


2 




Mg XII 


1 470+0.003 

-0.002 


4±2 


0.7 ±0.3 


2.3 


1.4697 ±0.0015 


3 (fixed) 


0.5 ± 0.3 


1.7 




Si XIV 


2.003 ± 0.002 


3.5 ± 1.3 


0.6 ± 0.2 


3 


2.002 ± 0.002 


3.5 (fixed) 


0.5 ± 0.2 


2.5 




Fei 


6.392 ± 0.005 


17 ±5 


2.7 ±0.5 


5.4 


6.383+«»« 


< 13 


1.5 ±0.4 


3.8 




Fexxvi 


6.958 ±0.013 


< 33 




2.6 


6.940 ± 0.015 


14 (fixed) 


0.8 ±0.5 


1.6 




VII (i) 


0.5686 ± 0.0003 


1.5 ±0.4 


28 ±3 


9.3 


0.5685 ± 0.0006 


< 2.1 


11 ±3 


3.7 




O VIII 


0.6521 ±0.0011 


3 (fixed) 


5.3 ± 1.1 


4.8 


0.653 ± 0.002 


3 (fixed) 


4.0 ± 1.0 


4 




Ne IX (i) 


0.917 ±0.003 


5+3 


3 0+°' 


6 


0.913 (fixed) 


7 (fixed) 


1.5 ±0.6 


2.5 


XMM 


Si XIV 


1.98 ±0.02 


68 ±24 


2.3 ±0.7 


3.3 


1.97 ±0.02 


60 (fixed) 


1.2 ±0.5 


2.4 




Fei 


6.398 ± 0.008 


47 ± 13 


3.3 ±0.3 


11 


6.40 ± 0.02 


108 ± 30 


2 6+° '' 
06^8:^ 

"•"-0.4 


8.7 




Fe XXV (i) 


6.66 ± 0.03 


20 (fixed) 


0.8 ±0.2 


4 


6.67 (fixed) 


20 (fixed) 


1.5 




Fe XXVI 


6.98 ± 0.02 


80 ±20 


1.7 ±0.3 


5.7 


6.97 ± 0.02 


65 ±30 


1.2 ±0.3 


4 








Id. 


E 
(keV) 


<T 

(eV) 


1(10-'*) 
(ph cm"^ s"') 


significance 







0.95-0. 1 .04 Phase-Interval 



HETGS 



XMM 



Ovn (i) 
O VIII 

Ne IX (i) 
Nex 

Mgxi(i) 
Mgxii 
Si XIV 
Fe I 

Fe XXVI 

Ovii(i) 
O VIII 
Ne IX (i) 
Si XIV 
Fe I 
Fe XXV (i) 
Fe XXVI 



0.569 (fixed) 
0.653 (fixed) 
0.921 ± 0.002 
1.020 ±0.002 
1 340+0-002 

■^"-0.003 

1.469 ±0.002 
2.005 ± 0.003 
6.378 ±0.011 

6 94+°"5 

"■^ -0.02 

5694+° ™*' 
0.653 (fixed) 
0.913 (fixed) 
2.003 (fixed) 
6.38 ± 0.03 
6.67 (fixed) 
6.93 ± 0.07 



1.1 (fixed) 

1.3 (fixed) 

2.4 (fixed) 

2 (fixed) 

3 (fixed) 
3 (fixed) 

3.6 (fixed) 
12 (fixed) 
25 (fixed) 

1.5 (fixed) 
3 (fixed) 
7 (fixed) 

60 (fixed) 

< 65 
20 (fixed) 
60 (fixed) 



< 140 
<26 
4±2 

1.5 ± 1.2 
1.0 ±0.5 
0.6 ± 0.4 
0.6 ±0.3 
1.2 ±0.6 
1.2 ±0.9 

20 ±5 
3±2 
1.5± 1.1 

< 1.2 
1.2 ±0.4 

< 0.6 
0.7 ± 0.4 



2 

1.3 

2 

1.5 

2 
2 

1.3 
4 

1.5 
1.4 



1.8 



Note — Best-fit values of the emission lines obtained using Gaussian profiles. Uncertainties are at the 68% c. 1. for each shown parameter. 



4. Discussion 

We analysed one XMM and two HETGS observations of XI 822- 
371, studying the averaged spectrum and the orbital-phase re- 
solved spectra. The best-fit value of the equivalent hydrogen col- 



umn associated with neutral interstellar matter, Nh ~ 0- 16 X 
10^^ cm"^, is similar to that obtained by llaria et al.l (1200 lb (i.e. 
0.11!°-^ X 1 0^^ cm-^) and to the value of ~ 0.16 x 10^^ ^m-^ 
obtained by 'J i et all (1201 Ih . who fitted the HETGS spectrum of 
XI 822-37 1 with a power-law component absorbed by local neu- 
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Table 7. Line equivalent widths along the orbital period 



0.04-0.25 0.25-0.50 0.50-0.75 0.75-0.95 0.95-1 04 

~~RGS Ovii (i) 24±"6 23T5 25T4 ]A±6 28^^^ 

RGS Ovm 4±3 6±2 6±2 7±4 <14 

HETGS Neix(i-r) 7 + 5 7±3 5±3 5±4 <14 

HETGS Nex 4±2 4±2 <4 <5 <7 

HETGS Mgxi(i) 2.8 ±1.1 4.3 ± 1.4 3 ±2 <3 <6 

HETGS Mgxii 2.6 ±1.1 1.8 ±1.0 <3 <3 <5 

HETGS Sixiv 3.3 ± 1.3 2.3 ± 1.3 2.2![;^ 2.1!];* <5 

HETGS Fei 50 ± 13 46 ± 12 36 ± 14 23 ± 14 < 46 

EPIC-pn Fei 54 ±11 62 ± 7 55 ± 9 50 ± 13 33 ± 20 

HETGS Fexxvi 2S+\l 20+\l < 36 < 30 < 67 

EPIC-pn Fexxvi 31 ± 10 24 ± 6 31!^- 25 ±11 < 47 

Note — The equivalent widths of the most prominent lines observed with HETGS, RGS, and EPIC-pn. The values are in units of eV, the errors 
are at 68% confidence level. 



low 0.8 keV. Future observations, e.g. with Astro-H, will allow 
one to asses the physical origin of this component. 



d 



0.50-0.75 
Orbital Phase 



Fig. 9. Best-fit values of the equivalent hydrogen column density 
associated with the local neutral matter (upper panel) and the 
fraction / of the emitting region covered by local neutral matter 
(lower panel). 



The 0.1-100 keV extrapolated unabsorbed fluxes of the 
Comptonised component obtained from the XMM and Chandra 
averaged spectra are 9.2 x 10"'" and 1.3 x lO"** erg cm"^ s"', 
respectively. To check the consistency of the measured fluxes in 
the two observations, we inspected at the long-term light curve 
provided by the All-Sky Monitor (ASM) on board the RXTE 
satellite. During the XMM and Chandra observations the ASM 
averaged count rate was 1.3 + 0.2, and 1.7 + 0.2 c/s, respec- 
tively. The ratio between the ASM count rate during the Chandra 
observation, ASMc/,, and the XMM observation, ASM^mm, is 
1.3 + 0.3, and the ratio of the 0.1-100 keV extrapolated unab- 
sorbed Chandra and XMM flux is 1.4, compatible with the ra- 
tio ASMca/ASMxmm- The proportionaUty of the ASM count 
rates to the unabsorbed fluxes suggests that the spectral shape of 
X1822-371 does not change significantly with time. Adopting a 
distance to the source of 2.5 kpc (Mason & Cordova 1982, here- 
after we will use this distance), the extrapolated observed lumi- 
nosity in the 0.1-100 keV energy range is 7 x 10^^ erg s"' and 
10^^^ erg s ' for the XMM and Chandra averaged spectra, respec- 
tively. 



tral matter The values of equivalent hydrogen column associated 
with local neutral matter, N//^^ , obtained by fitting the XMM and 
Chandra averaged spectra are compat ible with each ot her and 
consistent with the value reported by laria et all (1200 ih (about 
4.4 x 10^2 cm-2) and bv iJi et al.1 (1201 Ih (between 5.4 and 6.3 
xlO^^ cm"^). Finally, we found that the covered X-ray emitting 
surface is about 60-65%. 

The spectra show several emission lines associated with 
ionised elements and a fluorescence iron line, as already shown 
bylCottam et al. (2001) and Ji et al. (2011). 

We observed significant residuals in EPIC-pn data below 0.8 
keV. To remove them we added a black-body component with 
a temperature of 0.06 keV and a corresponding luminosity of 
(6.3 + 1.3) X 10^4 ej-g/s, a factor of one hundred smaller that the 
luminosity of the Comptonised component. We will not discuss 
this component in more detail because we cannot distinguish be- 
tween a real effect of an instrumental artefact of the Epic-pn be- 



In the following we propose a possible interpretation of the 
features observed in the X-ray spectrum of XI 822-371. We pro- 
pose a scenario in which 1) the modulation of the light curve 
is mainly caused by the occultation by the outer accretion disc 
(section l4Tb : 2) only 1% of the luminosity produced in the in- 
nermost region of XI 822-371 is scattered by an extended opti- 
cally thin corona into the line of sight (section l4!2l i: 3) the pho- 
toionised lines of Ovm, Neix, Nex, Mgxi, Mgxii, and Sixiv 
are produced in a compact region at the outer radius of the ac- 
cretion disc (section I4.4l i: 4) the Fe xxvi line is produced at a 
distance from the neutron star less than 3.7 x 10'" cm (section 
14.4b ; 5) the Ovii and the fluorescence iron line are produced in 
the photoionised surface of the accretion disc at radii between 
2 X 10'° and 4 x 10'° cm (sections |4i4] and |431i; 6) the neutral 
matter that partially covers the scattered emission is placed at a 
distance from the neutron star comparable to, or larger than, the 
outer accretion disc radius (section l43T l. 
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4.1. The presence of an opaque shield 

Our best-fit model consists of an optically thick Comptonised 
spectrum absorbed by neutral interstellar matter that is par- 
tially absorbed and Thomson-scattered by local neutral matter. 
However, even taking into account the Thomson scattering, this 
model does not explain the different values of the unabsorbed 
flux for the five selected phase-intervals. The 0.4-12 keV un- 
absorbed flux obtained from the XMM data changes along the 
orbital period from (7.22 + 0.03) x 10""^ erg cm"^ s"' in the 
phase-interval 0.25-0.50 down to (4.05 + 0.06) x l O''" erg cm'^ 
s~' du ring the eclipse (see Fig. [8] bottom panel). iBurderi et alj 
(120101) gave strong evidence that the neutron star in X1822-371 
accretes at the Eddington limit while the companion star has 
an overall mass outflow between 3.5 and 7.5 times the max- 
imum Eddington accretion rate riiE for a 1.4 M© neutron star 
Consequently, we should expect that the unabsorbed flux of the 
system is constant at the Eddington limit along the orbital period. 
The observed changes of the flux can be explained by the pres- 
ence of an opaque shield whose height varies along the orbital 
phase. The presence of an opaque shield i n XI 822-371 was al- 
ready discussed by several authors (see e.g. White & Holt '1982"; 
Mason & Cordova 1982; HelHer & Mason 1989; Bayless et al. 
20101) . who suggested that the outer accretion disc of X1822- 
371 is geometrically thick and its height depends on the az- 
imuthal angle. The 0.4-12 keV unabsorbed flux for the phase 
intervals 0.04-0.25, 0.50-0.75, 0.75-0.95, and 0.95-1.04 is 83%, 
91%, 76%, and 56%, respectively, as a percentage compared to 
that during the phase interval 0.25-0.50. The highest flux reduc- 
tion during the phase-interval 0.95-1.04 is due to the occultation 
caused by the companion star. 

On the other hand, the equivalent hydrogen column of the 
local neutral matter and the corresponding covering fraction / 
are correlated with the light-curve modulation (see Figs. |2] and 
|9|l in the sense that their best-fit values increase when the count 
rate increases and reach the minimum when the count rate is 
minimum (out of the eclipse). This suggests that we cannot as- 
cribe the observed light-curve modulation to the local neutral 
matter because if that were the case, we should expect a rise of 
the equivalent hydrogen column of the local matter and/or of the 
covering fraction / when the count rate drops at phases 0.75- 
0.95. 

Because we identify the opaque shield with the outer accre- 
tion disc, we expect that it is the highest when the unabsorbed 
flux is minimum and, conversely, the lowest when the unab- 
sorbed flux is maximum; in other words, the height of the outer 
disc is minimum at phases 0.25-0.50 and maximum at phases 
0.75-0.95. In agreement with lCottam et al.l (1200 1) and other au- 
thors we propose that the changes of the height at the outer rim is 
caused by the impact with the outer disc of matter coming from 
the companion star. At the impact point we expect the largest 
height of the outer radius of the accretion disc. 

4.2. The continuum emission and presence of an optically 
thin ADC 

Althou gh we find the same spectral shape reported bv llaria et al.l 
(1200 Ih . in this work we will give a diff'erent inte rpretation of 
its orig i n in light of th e recent result s reported by iBurderi et al.l 
JIqIO) . iBavless et al.l (l20Toh . and llaria et al.1 (l2oT l]), which 
showed a large derivative of the orbital period of X1822- 
371. Burderi et al. (2010) showed that the large orbital-period 
derivative of XI 822-371 indicates that the system accretes at 
the Eddington limit. Another indication that the luminosity 



of XI 822-371 is at the Eddington Umit is given by the ratio 
Lx/L„i„, where Lx and L„,„ are the X-ray and optical luminos- 
ity. fHeflier & Mason (1989) showed that the ratio Lx/Lo,,, for 
XI 822-371 is ~ 20, a factor 50 smaller than the typical value 
of 1000 for the other LMXBs. This indicates that the intrinsic 
X-ray lu minosity is underestimated b y at least a factor of 50. 
Finally, J onker & van der KiisI (1200 Ih ' showed that for a lumi- 
nosity of X1822-371 of 10^^ erg s~' the magnetic field intensity 
of the neutron star assumes an unlikely value of 8 x 10'^ G, while 
for a luminosity of the source at the Eddington limit it assumes 
a more conceivable value of 8 x 10'° G. These three indepen- 
dent arguments lead us to assume that the intrinsic luminosity of 
XI 822-371 is close to the Eddington limit; hence the first point 
that we discuss is why the observed lurainosity of XI 822-371 
is lO^ ^ *" erg s'' (see e.g.'Heflier & Mas onlll989l:lHeinz & NowakI 
l200ll:lParmaretaL 2000; laria et al. 2(3011 and this work). 

Considering that the outer disc is thick and that the inclina- 
tion angle of the system is 82.5° + 1.5°, we do not observe the 
direct emission from the innermost region and invoke the pres- 
ence of an extended corona above the disc that scatters 1% of 
the Eddington X-ray luminosity along the line o f sight. A sim- 
ilar scenario was discussed by McClintoc k et al.l (Il982h for 4U 
2129-1-47 and by Heflier & Mason (1989) for this source. This 
scenario is different from that proposed by laria et al.l (1200 ll) and 
other authors, according to which the continuum emission origi- 
nates in an extended optically thick corona that has a luminosity 
of 10^^ erg s '. Indeed, if an extended optically thick corona is 
present, a large part of the Eddington luminosity should be repro- 
cessed there and re-emitted, and the observed luminosity would 
be almost lO^*** erg s"', unless one considers the unlikely sce- 
nario in which the emission from the inner region of XI 822-371 
is extremely beamed. 

The photons produced in the innermost region near the NS 
are Comptonised in an optically thick plasma (t ^ 20) with 
a temperature of 3 keV, producing the main spectral compo- 
nent observed in the X-ray spectrum. The Comptonising cloud is 
probably a compact region near the NS. XI 822-371 has a pecu- 
liar geometry for two reasons: the system is observed at high in- 
clination angle and the height of the outer accretion disc is large. 
Combining these two eff'ects, we infer that we do not observe the 
direct emission Lq from the central region, but only the emission 
scattered by an extended corona. The scattered flux, considering 
Tc « 1, is Lq [(1 - e-^^)/i4nD^)\ ^ LoTc/iAnD-), where D is 
the distance to the source. Therefore the scattered flux is Tc times 
the emission produced near the NS. Assuming that tc ~ 0.01 we 
match the expected luminosity at the Eddington limit and the ob- 
served luminosity of ~ 10^^ erg/s. 

iJonker et al.l (l2003l) found that the mass function of X1822- 
371 is (2.03 +0.03) xIO^^Mq. Assuming a typical NS mass of 1.4 
Mq, we derive a companion star mass of 0.408 Mq. Considering 
that the orbital period of X1822-371 is 5.57 hours, we derive 
that the orbital separation is a = 1.34 x 10" cm using the third 
Kepler law. The Roche lobe of the NS and the companion star 
aieRi, = 9.8x10"' cm and/?/,, = 3.7xl0'° cm, respectively. For 
a mass ratio q = M2/M1 ^ 0.3, the tidal truncation radius, that 
is the maximum possibl e value of the accretion disc radius R^, i s 
0.43fl = 5.7 X 10'° cm (iFrank et al.ll2"ool iBavless et al.ll2010l) . 
We adopt Rd = 5.7 x 10'" cm as outer accretion disc radius. 
The estimated orbital parameters allow us to infer the size of 
the optically thin corona, making the reasonable assumption that 
the Roche lobe radius of the companion star is equal to the 
radius of the companion star R2. 
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We can estimate the size of the ADC taking into account 
the different count rates observed at mid-ecHpse time and out of 
the eclipse. Since the spectral shape does not change along the 
orbit, we impose that the change in count rates, in and out the 
eclipse, is mainly due to the volume of the visible optically thin 
ADC occulted by the companion star. Looking at the HETGS 
folded light curves of X1822-371 in Fig.|2l we observe that the 
count rate is Ced ~ 3 cts s"' at the mid-eclipse time while it 
changes between 4 and 6.6 cts s"' out of the eclipse. The maxi- 
mum observed count rate C„,„ indicates the minimum shielding 
and allows us to give a more accurate estimate of the ADC size. 
We also note that the values of C^ci and C„,„ in the EPlC-pn light 
curves are 37 and 84 cts s"', respectively; we obtain consistent 
results using EPIC-pn or HETGS count rates. Assuming a spher- 
ical volume of the ADC, the scattering volume is proportional to 
Cout out of the eclipse, while the scattering volume of the ADC 
subtracted by the volume occulted by the companion star is pro- 
portional to Ceci during the eclipse. That is 



3 ADC 



Cout 



2.2. 



We find that Radc ~ 6 x 10'" cm similar to the accretion disc 
radius R4 ~ 5.7 x 10'° cm. In the following we assume that 
Radc - Rd ~ 5.7 x 10'*^ cm. Considering the ADC as an ex- 
tended optically thin cloud with tc ~ 0.01, we can estimate the 
averaged electron density «g of the ADC using tq - (TTrieRhDC 
where crj- is the Thomson cross section and /?adc is the radius 
of the accretion disc corona; we obtain that the average electron 
density in ADC is ~ 3 x 10" cm"-'. 

The detection of the NS s pin period of 0.59 s in XI 822-371 
(iJonker & van der Klisll200l1) also supports a scenario in which 
an extended optically thin corona is present. Indeed, if we as- 
sume that the extended corona is optically thick and the observed 
Comptonised spectrum comes from the extended corona, then 
the extended corona should have an optical depth of r = 20, 
as observed in our fits. The unscattered photons fraction, F - 
Foe^'^, should be extremely small because of r ^ 20. All pho- 
tons in the corona should be scattered with a number of steps to 
escape of N = max(r, t^) ^ (for r » 1) and their travelled 
space should be NA = t^A, where A is the mean free path of the 
photons. For an optically thick ADC with Radc - 5.7 x 10'° cm 
and T ^ 20 we find n^, ^ 5 x 10'"* cm"-^, /I ^ 3 x 10^ cm, the 
average delay of the escaping photons - t^A/c ^ 40 s, and, 
finally, the spread in the arrival times of the coherent photons 
Af - tA/c ^ 2 s. Note that Af is almost three times the observed 
pulse period of XI 822-371 (that is 0.59 s), suggesting that a sce- 
nario with an extended optically thick corona in XI 822-371 is 
not viable because the ADC should wash out coherent pulsa- 
tions and prevent the observation of the 0.59 s pulse period. On 
the other hand, our scenario predicts an extended optically thin 
ADC with an electron density tig ~ 3x10" cm"^ and tc - 0.01 ; 
the mean free path of the photons in the ADC is /I ~ 5x10'^ 
cm, a factor 40 larger than the orbital separation of the system, 
in agreement with the fact that most photons are not scattered in 
the corona. We expect that only 1% of the photons are scattered 
in ADC, the escaping photons are scattered once at most and 
the average delay of the escaping photons is r/l/c ^ 2 s, similar 
to the optically thick case. However, because of the absence of 
self-occultation in ADC, the spread in the arrival times of co- 
herent photons depends on the travelled length I in ADC. In the 
worst case, for photons scattered in the outer region of the ADC, 
I = Radc then At„,ax - Radc I c ~ 2 s, while for photons scat- 
tered in inner regions of the ADC, for example at /? = 10'" cm. 



we find Af - 0.3 s, which allows us to observe the 0.59 s pulse 
period, although with a reduced rms amplitude. We note that the 
1 % scattered emission coming from the innermost region is ad- 
ditionally shielded by the outer disc, as discussed in the previous 
section. Finally, the different fluxes measured during the XMM 
and Chandra observation can be ascribed to a slightly different 
optical depth of the scattering corona or shielding from the outer 
disc rather than to an intrinsic variation of the sour ce luminosity 
that constantly emits at the Eddington limit (see iBurderi et al.l 
(2010) ). 

4.3. Presence of neutral matter around the binary system 

We detected the presence of local neutral matter that partially 
covers the X-ray scattering ADC. We found that the local neutral 
matter has an equivalent hydrogen column of around 5 x 10^^ 
cm"^ and a covering fraction of 60-65%. 

We propose that a large part of the outflowing matter is re- 
sponsible for the partial occultation of the extended X-ray scat- 
tering ADC, and suggest that the accreting matter above the 
Eddington threshold is radiatively expelled from the system at 
the inn er Lagrangian point. Indeed, as estimated bv lBurderi et al.l 
(20 Id), only the 20% of the mass transferred from the compan- 
ion star accretes onto the neutron star of 1.4 Mq. We report be- 
low our estimation on the equivalent hydrogen column associ- 
ated with the local neutral matter considering the mass continu- 
ity equation: 



{l-/3)\M2\=A7:r^^p(r)v(r). 



(1) 



In eq. [T] /? is the fraction of matter from the companion star that 
accretes onto the compact object, M2 is the rate of mass transfer 
from the companion star, ^ is a parameter that takes into account 
a non-spherical distribution of matter, p{r) and v(r) are the den- 
sity and speed of the matter at a given radius. The speed of the 
matter will be faster than the escape speed Vesc('"); and hence 
v(r) = ?7V'esc('") with 77 a parameter larger than 1. 

We know that the escape speed can be roughly written for a 
distance r » a (with a orbital separation of the binary system) 



G{Mi +M2) 



1/2 



(2) 



where Mi and M2 are the masses of the neutron star and com- 
panion star, respectively, and G is the gravitational constant. 
Combining eq.|2]with eq.[T] we obtain 

p(r) = 4.7x10^ (l-y8)r'?7"' mi, (mi+mzr^'^r-^^^ g cm-\(3) 

where is the mass transfer rate from the companion star in 
units of Eddington mass accretion rate and mi and m2 are the 
masses of the neutron star and companion star in units of solar 
masses. We obtain the particle density n from eq. [3] by divid- 
ing the density p by the proton mass and taking into account an 
appropriate fully ionised cosmic mixture of gases (fi = 0.615): 

n(r) 4.5 X 10^^ (1 -/3)r''7"' m2, (nii +m2y^'^r-^'^ cm"l(4) 

From the third Kepler law, the orbital separation a is 

a = 3.5 X 10'°(mi + otz)'^^^^^ cm, 

where Pi, is the orbital period in units of hours. We can rewrite 
eq.|4]as 

/ r\-3/2 

n(r) 6.9x10" (l-y6)r'??"' (mi -htoz)"'^',;' [-] cm-^(5) 
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iBurderi et alJ (l20ld) reported a value of ^ 0.18 and - 
5.5 for a neutron star mass of 1.4 Mq. The orbital period of 
X1822-371 is P/, ^ 5.57, therefore we derive 

n(a) ^ 3.1 X 10"(^77r' cm^l 

Supposing a constant particle density along the line of sight, 
we can determine the equivalent hydrogen column associated 
with the neutral matter using - n(a) x a, where a ^ 1.3 x 10" 
cmforX1822-371,thatis 

A? ~ 4 X 10^^(f77)"' cm"^ 

This value is similar to the equivalent hydrogen column of the 
cold matter that partially occults the central region 5 x 10^^ 
cm"^) under the hypothesis that the product (rj is close the unity. 
This agrees with our scenario, as ( takes into account the non- 
spherical distribution of matter and assumes values lower than 
1, while 77 is a parameter that indicates the speed of the matter 
in units of the escape speed, the lowest possible speed under the 
hypothesis that the matter leaves the binary system. The pres- 
ence of matter in the outer region of the disc or i n a region sur- 
roundi ng th e enth'e system was also proposed by iHarlaftis et al.l 
(Il997h and ICowlev et aP (l2003h . observing that the Ha Une is 
formed in a region not strongly eclipsed and that the Ho- veloc- 
ity variati£nJsj[eQ;_sm results were confirmed recently 
bv lPeris & Dil Vrtilekl (l2012h 

4.4. The ionised matter 

We detected ionised emission lines associated with He-like and 
H-like ions that are Ovii, Oviii, Neix, Nex, Mgxi, Mgxii, 
Si xm. Si XIV, and Fe xxvi. The Ne ix line is identified with the 
intercombination line during all orbital phases except for that 
during the eclipse (phase interval 0.95-1.04) when it is identified 
with the resonance line. We found that all the observed emission 
lines are clearly detected during the orbital phases between 0.04 
and 0.75. At the phase interval 0.75-0.95 we clearly detected the 
Ovii and Oviii in the RGS spectrum, of the neutral/partially- 
ionised iron line in the HETGS and EPlC-pn spectra, and, fi- 
nally, of the Fe xxvi line in the Epic-pn spectrum; the other lines 
have a significance lower than 3cr that does not allow us to state 
with certainty whether they are present at those orbital phases 
even though the small associated errors with the line energies 
(which are all compatible with the expected rest-frame energies) 
suggest their presence. At the eclipse (phase-interval 0.95-1.04) 
only the O vn line and the neutral/partially-ionised iron line are 
clearly identified in the RGS and EPlC-pn spectrum, respec- 
tively. 

The same em ission lines were already observed by 
ICottam et al.l (1200 l l). who analysed a shorter Chandra observa- 
tion. However, because of the poorer statistics, the Unes were 
observed only during the phases 0-0.5. ICottam et all ('2001) sug- 
gested that the lines originate from the inner face of the bulge 
at the outer radius of the accretion disc, where the stream of ac- 
creting matter impacts the outer region of the disc. During the 
phases 0-0.5 we directly observed the inner face of the bulge, 
whereas in the phase-interval 0.5-1 the lines are not observed 
because the bulge is facing the opposite direction with respect to 
the observer. 

We report the equivalent widths of the most prominent lines 
at different phases in Table|2] The equivalent widths of the emis- 
sion lines associated with O viii, Ne ix, Ne x, Mg xi, Mg xii, and 
Si XIV do not show significant changes at the phase intervals 



where they are clearly detected, implying that the line and con- 
tinuum fluxes change in the same way. The equivalent width of 
the O VII seems to show a reduction at phase 0.75-0.95, and the 
equivalent width of the Fe xxvi line does not vary along the or- 
bit, although we have no clear detection of this line during the 
eclipse and more investigations at this phase interval are needed; 
however, this line cannot originate in the same region of the 
bulge because it has largest associated ionisation parameter with 
respect the other lines. 

The He-like triplets of O vii, Ne ix, and Mg xi show that the 
intercombination lines are dominant with respect to the forbid- 
den and resonance lines at all phases. In the following we re- 
port our plasma diagnostics using the best-fitting values obtained 
from the averaged spectrum, because they provide the best con- 
straints. From the line intensities of the O vii triplet, we infer that 
the origin of the lines is a photoionised plasma with an electron 
density n,, higher than lO'^ ^^n^ (see the bottom right panel in 
Fig. 11 in iPorquet & Dubaij||2000 ). Our analysis indicates that 
the R parameter (ratio of the forbidden over intercombination 
line intensity) is 7? < 0.2 for O vii, Ne ix and Mg xi, suggesting a 
plasma density > 10'^ cm~^ for O vn, > lO'^ fo r Ne ix, and 
> 10'"^ cm"^ for Mgxi (see lPorquet & Dubaiill2000h . Assuming 
that the Ne ix and Mg xi are produced in the bulge, we derive that 
the bulge has an electron density > lO'"* cm"^. 

iKallman & Bautista (2001) studied the properties of a pho- 
toionised plasma at high density (lO"" cm"^) and derived that 
the ionisation parameter ^ is between 100 and 400 for Neix, 
Nex, Mgxi, Mgxii, and Mgxii and Sixiv, respectively, and 
^ ^ 1000 for Fexxvi. We can estimate the electron density 
of the photoionised plasma using the relation ^ = LxlifieP') 
(see Krolik et al. 1981), where Lx is the X-ray luminosity of the 
source, ^ the ionisation parameter, the electron density of the 
matter, and r the distance of the photoionised matter from the 
central source. Adopting an X-ray luminosity of 2 x 10""^ erg/s, 
and placing the bulge at 6 x 10"^ cm (the outer radius of the ac- 
cretion disc), we derive an electron density ^ 5 x lO''* cm"^ 
that is compatible with the value of > lO'"* cm"-' estimated 
above from the study of the triplets. 

To estimate the temperature of the emitting plasma we use 
the G parameter inferred in section [372] which is 4.3;^j 5, 4.4^j 3, 
and 3.7;^j^ for O vii, Ne ix a nd Mg xi, respectiv ely. According to 
the calculations outlined in iPorquet & Dubaul ([2000), we deter- 
mine that the temperature of the plasma should be < 5 x 10^ K, 
< 8 X 10^ K, and < 10^ K for O vii, Ne ix and Mg xi, respec- 
tively. ICottam et al.l (1200 Ih detected a RRC feature associated 
with Ne IX with a temperature of 1 3 ± 7 e V corresponding to a 
plasma temperature of (1.5±0.8)xlO^ K. We looked for a similar 
feature in our HETGS observation and found only an upper limit 
to the plasma temperature of 14 eV corresponding to 1 .6 x 10^ K. 
We conclude that the region where the Ne ix, Ne x, Mg xi, Mg xn, 
and Mg xn and Si xiv lines originate has an electron density of 
Hf, - 5 X 10'"* cm"^ and a tem perature of ~ 1.6 x 10^ K. 

Recently, Ji et al.l (1201 ih . analysing the same Chandra data 
sets as used for this work, measured a Doppler-shift of several 
lines. Iji et al.l(l201 ih divided the folded light curve into 40 phase- 
intervals (separated by 0.025 in phase) and adopted a phase win- 
dow of 0.16 width in phase, sliding along the 40 intervals from 
which they extracted the spectra. In this paper we follow a diff'er- 
ent approach, obtaining information about the lines from the five 
independent phase-intervals discussed above. Below we present 
the Doppler-shift velocities of the O vn line observed with the 
RGS instrument, and of the Ne ix, Ne x, Mg xi, Mg xn, and Si xiv 
lines observed with the HETGS instruments. We estimate the 
Doppler-shift velocities of the O vii, Ne ix, Ne x, Mg xi, Mg xn. 
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Fig. 10. Best-fit parameters of the Doppler modulation associated with the O vii, Ne ix, Ne x, Mg xi, Mg xii, and Si xiv emission 
lines adopting a sinusoidal function as discussed in the text. Uncertainties are at 68% c. 1. 

Table 8. Best-fit parameters of the Doppler-shift modulation of the most prominent emission lines. 



Line 


Rest-frame energy 




K 


<Po 


x'(d-o.f) 


Prob. chance, improv. 




(keV) 


km/s 


km/s 






(%) 


Ovii 


0.5687 


-44 ± 87 


77 ± 121 


0.20 (fixed) 


1.8(3) 


52.6 


Ne IX (i-r) 


0.9149(1); 0.9220(r) 


-11 ±97 


385 ± 165 


0.22 ± 0.08 


0.012(2) 


99.9 


Nex 


1.0218 


-95 ± 57 


330 ± 77 


0.32 ± 0.05 


1.51(2) 


94.1 


Mgxi(i) 


1.3433 


-60+152 


306 ± 122 


0.24 ±0.10 


0.89(2) 


88.7 


Mgxii 


1.4723 


-303 ± 143 


340 ± 153 


0.17 ±0.08 


0.29(2) 


94.7 


Si XIV 


2.0055 


-217 ± 113 


383 ± 150 


0.14 ±0.06 


1.82(2) 


78.5 


Fe I 


6.400 


-396 ± 104 


323 ± 155 


0.07 ± 0.07 


1.74(2) 


71.9 


Fe XXVI 


6.966 


-608 ±210 


235 ± 290 


0.83 (fixed) 


2.07(3) 


60.0 



Note — The parameters v„,„, K, and 0o are defined in the text. The probability of chance improvement is obtained with respect to a constant 
function. The errors are at 68% c.l. 



and Si xiv lines for each phase-interval using the best-fit values 
of the line energies obtained from RGS data for the O vii line 
and from HETGS data for the other lines, using the best-fit val- 
ues of the line energies shown in Table |6] This analysis was not 
possible for the O viii line because of larger uncertainties. 

To estimate the Doppler-shift modulations we fitted the 
Doppler-shift velocities, Vj/„y,, obtained for each phase-interval 
using the relation 

VMft = V,,.,, + K sm(2n((/) - (f>o)/P), (6) 

where the free parameters are Vjys,, the systemic velocity in units 
of km/s, K, the semi-amplitude velocity in units of km/s, and (f>o, 
the orbital phase at which the sinusoidal function is null. The 
orbital period P in units of phase is kept fixed to 1 . The best-fit 
parameters for each line are shown in Table [8] 



Initially, we focus our attention on the sinusoidal modulation 
of the Ne ix, Ne x, Mg xi, Mg xii, and Si xiv lines; the best-fit val- 
ues of Vsysr, K, and 0o are similar for all these lines, suggesting 
that these five lines are probably produced in the same region. 
We plot the obtained best-fit parameters of the sinusoidal mod- 
ulation in Fig. [TT] (left, middle, and right panel, respectively). 
Under the assumption that the emitting region is the same for 
the five lines, we fit their values of Vj,.,,, K, and (po using a 
constant and finding self-consistent results. We obtain in this 
way the following average values: Vsyst = -107 ± 41 km/s and 
x\d.o.f.) = 4.5(4),/: = 339 ±53 km/s andx^id.o.f.) = 0.25(4), 
and, finally, 0o = 0.23 + 0.03 and x^id-o.f.) = 6.1(4). 

The systemic velocity, Vsxsi, is similar for the five lines with 
a value of v,,,,, - -107 +41 km/s (see Fig.[TTl left panel), which 
is compatible with the systemic velocity of -106 km/s and -101 
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l'(dof).4.5(4) 
= -107 + 41 km/s 



I>(d.o.n=0.25(4) 
K=339±53 km/s 



Z"(d.o.0=6.1(4) 
phase-0.23±0.03 



NeX MgXli MgXll Si XIV 



Mg XI i Mg XII Si XIV 



Mg XI i Mg XII 



Fig. 11. Best-fit values of the parameters v,^.,, (left panel), K (middle), and (f>o (right panel) for the Neix, Nex, Mgxi, Mgxii, and 
Si XIV lines. Uncertainties are at 68% c. 1. 



km/s obtained by ICowlev et all (l2003h for the He n /14686 and 
Ha lines observed in XI 822-371. 

The values of 00 obtained for the five lines are similar, and 
fitting them with a constant, we obtain (pQ - 0.23 + 0.03 and 
a xlgjid-o.f.) - 1.5(4), which is acceptable. We find that the 
emitting region has null radial velocity at phases 0.23 + 0.03 and 
0.73+0.03 which is when the radial velocity changes from nega- 
tive to positive and from positive to negative values, respectively. 
The detection of a Doppler modulation implies that the Ne ix, 
Ne X, Mg XI, Mg XII, and Si xiv lines are produced in a compact 
region in the bulg e at the outer accretion disc, as suggested by 
ICottam et al.l (2001). Moreover, the phases between 0.6 and 0.9 
correspond to the minimum of the modulation out of the eclipse 
in the light curve (see Fig. |2]i, which is when the outer accretion 
disc is the tallest and, consequently, when the scattered emission 
is mostly occulted by the outer disc. Our results suggest that the 
outer region of the disc with the largest height and the bulge is 
the same region as the one that covers a phase interval longer 
than 0.3; while the region where the Ne ix, Ne x, Mg xi, Mg xii, 
and Si xiv lines are produced might be (a part of) the bulge illu- 
minated by the central source. 

We show the inferred geometry of X1822-371 in Fig. [T2j 
the region in the bulge where the lines originate (indicated with 
a white circle) is not to scale because its size is uncertain. 
Adopting the estimate of the solid angle subtended by this region 
inferred by Cottam et al. (2001) and the outer radius of the disc 
(6 X 10'° cm) estimated above, we find an illuminated surface 
area of ~ 8.3 x 10^" cm^. Assuming, arbitrarily, that the height h 
and the subtended arc I of the illuminated region have the same 
size, we find h = I ~ 3 x 10'" cm; the subtended angle in this 
case is ~ 28° corresponding to a phase interval of 0.07 that is 
shorter than that subtended by the whole bulge, which is 0.3. On 
the other hand, assuming that the whole inner face of the bulge 
is illuminated by the central source, we estimate that the vertical 
scale height is /i ~ 7.3 x lO** cm. Using detailed modelling of th e 
X-ray, UV, optical , and I R lig ht curves. White & Holtl ' (Il982h . 
iMason & Cordoval (|1982|) . and iHelher & MasonI (Il989l) all in- 
ferred a phase-dependent vertical structure along the outer ra- 
dius of the disc that is dominated by the bulge with a maximum 
vertical height between 0.6 and 1.6 xlO'" cm depending on the 
model. This is consistent with the vertical height that we esti- 
mate assuming that the whole inner face of the bulge emits the 
observed lines, while a discrepancy of a factor 2 is obtained for 
the first case (h ^ I). 

At phase 0.73 the outer part of the illuminated compact re- 
gion is facing the observer and changes its radial velocity from 



positive to negative values. At the same phase the outer region of 
the disc is the tallest and occults a large part of the direct emis- 
sion coming from the inner disc surface. The black-coloured arc 
at the outer radius of the disc shown in Fig. [T2]indicates where 
the outer disc is the tallest, that is the whole size of the bulge. 

The semi-amplitude velocities are very similar for the five 
lines, we find K - 339+53 km/s (see Fig.[TTl middle panel) with 
a x^id.o.f.) = 0.25(4). The maximum of the blueshift velocities 
(~ 360 km/s) is observed at phases betwe en 0.45 and . 5 1 an d 
is compatible with the value inferred by ICottam et al.l (1200 lb . 
who deduced a blueshift velocity of ~ 360 km/s at phase 0.50 
and suggested that we are observing the line emissions from the 
bulge that was e xpected to form in the shock-heated, colliding 
material (see e.g. lLivio et al. I IT986h . Our results suggest that the 
lines associated with the Ne ix, Ne x, Mg xi, Mg xii, and Si xiv 
transitions are produced in the bulge illuminated by the central 
source. 

The O VII line does not show a significant sinusoidal mod- 
ulation. We obtain that K - 11 + 121 km/s and the probabil- 
ity of chance improvement with respect to a constant is only 
52.6%. The flux and the equivalent width of this line drop dur- 
ing the phase interval 0.75-0.95 (see Tables|6]and|7]i. These argu- 
ments suggest that the O vii line is produced in a different region. 
Because the equivalent width of the O vii line drops at phases be- 
tween 0.75 and 0.95, we argue that at those phases its emitting 
region is occulted by the bulge and the outer disc, implying that 
the O VII line is produced in an inner region with respect to the 
other ionised lines. We can speculate that the Ovn line is pro- 
duced on the accretion disc surface illuminated by the central 
source. The half width at half maximum (HWHM) of the line 
Ovii corresponds to a velocity of about 680 + 130 km/s, and, 
consequently, to an accretion disc radius of (4.0;^j j) x 10'° cm. 
We estimate the electron density of the photoionised plasma as- 
suming a luminosity at the Eddington limit (2 x 10-"^ erg s"') 
using the relation ^ - LxlifieP') and obtain a plasma density of 
about 1 .3 X 10'^ cm"-', which agrees with our estimation from the 
study of the O vii triplet that «e is larger than lO'^ cm"^. The ra- 
dius where the O vii line originates is comparable with the com- 
panion star radius (that is 3.7 x lO"* cm), which explains why 
the equivalent width of the O vii line drops at phases between 
0.75 and 0.95 and increases again during the eclipse. Indeed, 
at phases 0.75-0.95, the bulge is interposed between the central 
source and the observer and its large azimuthal extension covers 
the region of the disc that emits the line, while during the eclipse 
this region is only marginally occulted by the companion star 
The O vn line is clearly detected in the RGS spectrum along the 
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Phase 



Phase 0.23 



Phase 0.73 





Observer 



Fig. 12. Illustration of the geometry of X1822-37. The figure is to scale except for the size of the region in the bulge where the 
Oviii, Neix, Nex, Mgxi, Mgxii, and Sixiv lines originate (white circle). The black-coloured arc at the outer radius of the disc 
indicates where the outer disc is the highest, which is the location of the bulge. The dashed and dotted circumferences on the disc 
indicate where the O vn and Fe i-Fe xv lines are produced, respectively. 



whole orbital period. Because we observe this line at phase in- 
terval 0.75-0.95, it is probable that we are observing both direct 
and scattered line photons of O vn. Indeed, if we had observed 
only the direct line photons, we should not have detected the 
O VII line at this phase interval. 

Furthermore, even if the emission lines associated with the 
Ne IX, Ne x, Mg xi, Mg xii, and Si xiv transitions are not clearly 
detected in the phase interval 0.75-0.95, we tentatively suggest 
that we are observing both direct and scattered emission from the 
illuminated region of the bulge; this scenario seems reasonable 
considering that the optically thin extended corona extends to 
the outer disc, implying that even when the outer face of the 
bulge is orientated towards the observer, we expect that part of 
the emission from the illuminated region is scattered into the line 
of sight by the extended optically thin corona. However, more 
investigations with observations with larger statistics are needed 
to confirm this scenario. 

We show in Fig. [12] the geometry of XI 822-371 at phase 
(during the eclipse) and at phase 0.73 (when the shielding by 
the bulge is maximum); during the eclipse the companion star 
occults the region of the disc where the O vn line was produced 
(indicated with a dashed circumference in Fig.[T2li less than the 
bulge occults the same region at phase 0.73. 

We were unable to estimate the Doppler motion associated 
with the O viii emission line because of the large errors associ- 
ated with its centroid energy. However, we suggest that this line 
is produced in the illuminated face of the bulge for two reasons: 
1) its equivalent width does not change along the orbital phase 
(see Tab. |7]i similarly to the equivalent width of the Ne ix, Ne x, 
Mgxi, Mgxii, and Sixiv emission lines; 2) the io nisation pa- 
rameter of the O VIII is near 100 (see lKallman & Ba utista 2001), 
similar to that obtained for the lines produced in the bulge (i.e. 
between 100 and 400). 

The Fe xxvi emission line does not show Doppler modula- 
tion. Initially, we fitted its Doppler-shift velocities with a con- 
stant and then with a sinusoidal function, obtaining a probability 
of chance improvement of only 60%. The corresponding equiv- 
alent width is constant at the phases 0.04-0.95 and, during the 



eclipse, it shows an upper limit compatible with the values ob- 
tained at the other phases (see Table|7]i. We have only a marginal 
detection of the Fexxvi em ission line during the eclipse. In 
agreement with Iji et al.l (1201 1) . we argue that this line may be 
produced at r < 3.7 x lO'" cm (i.e. the size of the Roche lobe 
radius of the companion star). 

4.5. The neutral/partially ionised iron line 

We find an emission line at 6.4 keV in HETGS and Epic-pn 
spectra, consistent with fluorescent emission from Fei-Fexv 
(Kall rnan et a l. 2004) with an ionisation parameter ^ < 10, sim- 
ilar to the one associated with O vii. The broadening of the line 
at 6.4 keV is 18 + 3 eV as measured by the HETGS (see Table 
|2|l. The width obtain ed from t he HE TGS data is compatible with 
the value inferred bv iJi et al.l (1201 Ih . 

We do not significantly detect Doppler modulation of this 
emission line (see Table |7]l. Its flux is minimum at phases be- 
tween 0.75-1.04, during the occultation by the bulge and the 
companion star. This can be explained by assuming that the flu- 
orescence iron line is produced at the surface of the accretion 
disc similarly to the O vii line, taking also into account that the 
ionisation parameter is comparable for the two lines. Using the 
HETGS best-fitting parameters and assuming that the broaden- 
ing of the neutral iron line is caused by a Keplerian motion of 
the plasma, we find that it should originate at a distance from the 
neutron star of (1.9^Q g) x 10'" cm, that is at a radius a factor of 
two smaller than that at which the O vn line is produced. Since 
the radius of the companion star is 3.7 x 10'" cm, the emitting 
region of the iron line is occulted both during the eclipse (phase 
0) and during the occultation by the opaque shield (phase 0.73). 
We show in Fig.[T2]the radius where the fluorescence iron line is 
produced with a dotted circumference. 



5. Conclusions 

We analysed one XMM-Newton observation and two Chandra 
observations of XI 822-371. We found that the continuum emis- 
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sion is well-fitted using a Comptonised component partially ab- 
sorbed by local neutral matter. In the following we summarise 
our results: 

- We inferred the presence of an extended optically thin 
corona. The Comptonised component, with electron temper- 
ature of 3 keV and optical depth of 20, is produced in a com- 
pact region near the neutron star In line with the work of 
iBurderi etaP (l2010h . we assume that the intrinsic luminos- 
ity is 2 X 10^^ erg/s (the Eddington luminosity for a neu- 
tron star of 1 .4 Mq) although the observed luminosity is lO-'^ 
erg/s. This discrepancy is caused by the high system inclina- 
tion and by the geometrically thick outer accretion disc that 
occults the innermost region. We observed only the emission 
that is scattered into the line of sight by an extended optically 
thin (t ^ 0.01) corona. 

- The partial covering component is due to neutral matter lo- 
cated at a radius larger than the outer disc radius. The frac- 
tion of the covered emitting surface is 60% along the whole 
orbit and the equivalent hydrogen column is between 4.5 
and 5 x 10^^ cm"^ . Since the mass transfer in XI 822-37 1 
IS not conservative ( iBurderi et alJl2010l:lBavless et alj|2010h . 
we suggest that the local neutral matter could be the matter 
outflowing from the companion star, which is swept away by 
the radiation pressure in the proposed scenario. 

- We observed the emission lines associated with Ne ix, Ne x, 
Mg XI, Mg xii. Si xm, and Si xiv. These lines show Doppler 
modulation, suggesting that they are produced in a compact 
region. In agreement with previous authors, we found that 
these lines are produced in the inner illuminated face of the 
bulge at the outer radius of the accretion disc distant 6 x 10'*' 
cm from the NS. The density of the photoionised emitting 
plasma is ~ 5 x lO''* cm"^. We detected the Oviii emis- 
sion line in the RGS spectrum at all phase intervals with 
a marginal detection during eclipse. The large errors on its 
centroid energy do not allow us to find a clear Doppler mod- 
ulation of this line. However, the behaviour of its equiva- 
lent width, similar to those of the Ne ix, Ne x, Mg xi, Mg xii. 
Si XIII, and Si xiv lines, and its similar ionisation parameter 
suggest that the O viii emission line is produced in the illu- 
minated bulge as well. 

- The Fe xxvi line was observed at all phases with a marginal 
detection during the eclipse. This line has a large ionisation 
parameter and does not show signific ant Doppler modula- 
tion. In agreement with Ji et al. I (I20TT1) . we suggest that this 
line is produced at a distance from the compact object less 
than 3.7 x 10'° cm. 

- We observed the O vii and a prominent fluorescence iron line 
that can be associated with Fei-Fexv. These lines do not 
show significant Doppler-modulation and are produced in 
the photoionised surface of the accretion disc, at 4.0;^j *j* x 
10'° cm and 1.9^q^ x 10'° cm, respectively. The O vii emis- 
sion line is partially occulted only at phases between 0.75- 
0.95 by the bulge while the fluorescence iron line is occulted 
also during the eclipse considering that the radius of the com- 
panion star is 3.7 x 10'" cm. 
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